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Abstract 

We present four-loop results for the gauge /3-function and the fermion mass 
anomalous dimension for a gauge theory with a general gauge group and a multiplet 
of fermions transforming according to an arbitrary representation, calculated using 
the dimensional reduction scheme. In the special case of a supersymmetric theory 
we confirm previous calculations of both the gauge /3-function and the gaugino mass 
/3-function. 
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1 Introduction 



In recent papers some of us presented calculations of the QCD /3-function, (3 S , and the 
fermion mass anomalous dimension (or mass /3-function), 7 m through three loops [1] 
and four loops [2] using the DRED (or DR) scheme, which is based on regularisation by 
dimensional reduction [3,4]. An interesting feature of these calculations is the dependence 
of /3 S on the evanescent couplings: e-scalar interactions that do not renormalise like the 
gauge coupling. At three loops j3 s depends on the e-scalar Yukawa coupling, and at four 
loops it also depends on the e-scalar quartic interaction. The first explicit calculations of 
the one loop corrections to this quartic interaction appeared in Ref. [5] (for a particular 
SU(2) model), and in Ref. [2] (for QCD). Here we generalise the calculation to SU(N), 
SO(N) and Sp(N). This involves some quite interesting and (relatively) little known 
group theory. We also similarly generalise the result from Ref. [2] for the e-scalar Yukawa 
coupling. 

Of course it is the SU(3) case described in the previous papers which is most obviously 
currently useful, but the general result is also of interest, for possible future applications 
to other symmetry groups, and if only as a further test of the validity of the DRED 
procedure. Our confidence in this is reinforced by once again comparing our results for 
the special case of supersymmetry (when there is a single fermion multiplet in the adjoint 
representation). The result for (3 S for a renormalisable M = 1 super symmetric theory was 
given through four loops in Ref. [6], the derivation being based on the completion of a 
construction of the coupling constant redefinition connecting the DRED scheme to the 
NSVZ scheme developed in Ref. [7] . Here we not only verify this result through four loops 
(in the special case of a theory with no superpotential) but also we verify the result for 
the gaugino /3-function through the same order. This is of interest because of course the 
gaugino mass breaks supersymmetry, and the issue of regularisation and renormalisation of 
softly-broken supersymmetric theories present additional subtleties. The exact formula for 
the gaugino /^-function (expressing it in terms of (3 S ), as first derived in Ref. [8] (inspired by 
an observation by Hisano and Shifman in Ref. [9]) relied heavily on the spurion formalism, 
as developed in particular by Yamada [10]; it is reassuring to find that the relationship 
between the two /3-functions indeed holds in an explicit DRED calculation. 

In section [2] we review the renormalisation procedure for a gauge theory using DRED; 
then in section [3] we describe the one loop renormalisation of the e-scalar self-interaction. 
We first give results for the SU(N) case, explaining how to reduce to the special cases 
N = 2 and N = 3. We then generalise to expressions valid for an arbitrary groups. 

In section H] we give the full four-loop results for f3 s and j m for the general case, and 
in section we reduce to the special case of supersymmetry for comparison with earlier 
results, as described above. Finally in the Appendices we explain some of the group 
theory involved in the calculations and give explicit results for SO(N) and Sp(N) for the 
one-loop e-scalar quartic interaction /3-functions. 
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2 Gauge theory with fermions 



Consider a non-abelian gauge theory with gauge fields and a multiplet of two- 
component fermions ijj A {x) transforming according to a representation R of the gauge 
group Q. 

The Lagrangian density (in terms of bare fields) is 

Lb = ~\G\ V - U d " W rf + CrPDfC? + i^ aA a^ a {D,) A B ^ (1) 

where 

G% = W - + gf bc WlW c v (2) 

and 

{D^) A b = S A B d, - ig{R a ) A B W; (3) 

and the usual covariant gauge fixing and ghost (C, C*) terms have been introduced. As 
usual cr M = (/, — <x) where cr are the Pauli matrices. 

For the case when the theory admits a gauge invariant fermion mass term we will have 
L B — > L B + L B , where 

T,™ = 

2 



Lb = \m AB r A € + cc. (4) 



Dimensional reduction amounts to imposing that all field variables depend only on 
a subset of the total number of space-time dimensions; in this case d out of 4 where 
d = 4 — 2e. We can then make the decomposition 

W;{xi) = {W?{xi),W a a {xi)} (5) 

where 

5\ = 5 j j = d and 5 aa = 2e. (6) 
It is then easy to show that 

Lb — L d B + U B (7) 

where 

L B = ~\G% ~ H&Wi? + C*d i D i C + h\xfD^ (8) 

and 

U B = i(AW CT ) 2 - g^R^W: ~ \g 2 f abc F de W b a W c a ,W*W e a ,. (9) 

Conventional dimensional regularisation (DREG) amounts to using Eq. ([HD and discarding 
Eq. ©• 

We would now like to rewrite Eq. and Eq. 1^ in terms of renormalised quantities. 
It is clear, however, from the dimensionally reduced form of the gauge transformations: 

6W? = diA a + gf abc W t b A c (10a) 
m a a = gf abc W b A c (10b) 
S^ A = ig(R a ) A B*p B A a (10c) 
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that each term in Eq. (Q is separately invariant under gauge transformations. The W a - 
fields behave exactly like scalar fields, and are hence known as e-scalars. There is therefore 
no reason to expect the ipif)W a vertex to renormalise in the same way as the ifiipWi 
vertex (except in the case of supersymmetric theories). In the case of the quartic e-scalar 
interaction it is evident that more than one such coupling is permitted by Eq. (llObj) . In 



other words, we cannot in general expect the f — f tensor structure present in Eq. (Q to 
be preserved under renormalisation. This is clear from the abelian case, where there is no 
quartic interaction in L e B but there is a divergent graph at one loop from a fermion loop. 

We are therefore led to consider the following expressions for renormalised quantities 
L d and L £ : 

L d = _ l z ww {diW ._ d . Wi) 2_j_ {d i Wi) 2 

_ z www gf ahc d i W a W u W ci - \Z m g 2 f abc f ade W b W-W di W ej 
+ Z cc d i C*d l C + Z ccw gf abc d l C a *W b C c 

+ ZWqtfdrt + ZWWg^K&TlW? (11) 

and 

re _ 1 ryee 
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Z ££ (diW a ) 2 + Z ££W gf abc diW^W bi W^ 



+ Z ££WW g 2 f abc f ade W b W^W di W^ - Z^ £ g e ^R a a a i\)W a a 
p 

- lJ2 Z r £X r H r bCdW aW^W b W d ,. (12) 



r=l 

In the case when we have a fermion mass term we would also have 



\Z m Z^mil) a ^ a + c.c. (13) 



Eq. (llip is the usual expression for the Lagrangian in terms of renormalised parameters. 
In Eq. ( fl2!) we have introduced a "Yukawa" coupling g e and a set of p quartic couplings A r . 
(Strictly speaking, Eq. (fl2l) should also have a mass term for the e-scalars; but since this 
mass term does not affect (3 S or 7 m we omit it here.) The number p is given by the number 
of independent rank four tensors H abcd which are non-vanishing when symmetrised with 
respect to (ab) and (cd) interchange. In the next section we discuss the quartic vertex 
and its renormalisation in more detail. 

1 Since e-scalars are present only on internal lines we could, in fact, choose the wave function renor- 
malisation of W a and W% to be the same; or, indeed, have no wave function renormalisation for W„ at all. 
The crucial thing is correct treatment of sub-graphs, which means recognition that vertices with e-scalars 
renormalise in a different way from their gauge counterparts. However, we choose to renormalise the 
e-scalar conventionally. 
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3 The e-scalar self coupling 



Let us discuss the structure of the quartic e-scalar couplings for an arbitrary gauge group. 
These interactions are invariant under the symmetry Q <g> 0(2e), where only the Q is 
gauged. The renormalisation properties of scalar theories with invariances of the type 
Gi <g> G 2 have been studied in considerable detail, for example 0(m) ® 0(n) in the theory 
of critical phenomena and U(m) <g> U(n) in the context of QCD. In these cases, however, 
the scalars transform as vector (fundamental) representations of the gauge group factors 
whereas for us they transform as adjoints. 

This raises an interesting group theory question: how many independent couplings 
are there for a given gauge group Q7 Evidently the question of how many independent 
tensors of the form K abcd there are is the question of how many times the singlet rep- 
resentation occurs in the reduction to irreducible representations of the product of four 
adjoint representations. (Neither this question, nor the obvious generalisation to n-tensors 
K ai '" a " has been much studied in the literature; an exception being the classic work of 
Cvitanovic [11], to which we will return presently). The set of tensors relevant to our 
problem is the subset of such tensors H abcd which is invariant with respect to (a, b) and 
(c, d) exchange, because of the 0(2e) invariance. 

If we have an irreducible basis of dimensionality 7(71) for the n-tensors of the form 

K^~ an , l<a< 7(71) (14) 
then a general n-tensor K ai "' an can be expressed in terms of the basis as 

K* 1 "*" =x f} K% 1 - an , (15) 
where xp are determined by the equation 

QlpXp = y a = K a - (16) 

and 

Q n a p = K 1 '" anK T" an - ( 17 ) 

Thus construction of the Q n -matrix permits reduction of an arbitrary n-tensor to the 
basis. 

3.1 The case Q = SU(N) 

The fundamental representation T a of the generators R a of SU(N) satisfies 



Tr (T a T b ) = V b , (18) 
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n 


N = 2 


N = 3 iV = 4 


iV = 5 


iV = 6 


2 


1 


1 1 


1 


1 


3 


1 


2 2 


2 


2 


4 


3 


8 9 


9 


9 


5 





Table 1: Basis dimensionality for n-tensors in SU(N). 



where b is a constant. For the rest of this section we will adopt the usual convention 
whereby 6 = 1. 

In table[T]we present some results for the dimensionality 7(71) for SU(N) as a function 
of N. It is interesting that Cvitanovic [11] remarks that a formula for the dimensionality 
of a basis (in general over-complete) is provided by the subfactorial (3{n) where 

/?(n)=n!(l-i + i + ...(-ir^). (19) 
1! 2! n! 

It appears that for sufficiently large N we have 7(71) = P(n). 



A natural choice for the basis for the case n = 4 when > 4 is given 




K 2 

The reduction of the basis to 7 = 8 in the case SU(3) is achieved via the relation [12,13] 

K 4 + K 5 + K 6 = l(K! + K 2 + K 3 ) (21) 

which is not valid for N > 4. The corresponding identity for general iV reduces a sym- 
metrised (A r +l)-tensor consisting of N—l (i-tensors; for an elegant derivation see Ref. [14]. 

2 An alternative way to define a basis which has the virtue of being immediately generalisable to 
any group [11] is in terms of traces of products of the generators in the defining representation, thus 
Tr (T Q T fc T c T' i ), Tr (T a T b ) Tr (T c T d ) etc. 



= S ab 5 cd 


K 4 


= d abe d cde 


K 7 


fjp-be jcde 


= 5 ac 5 bd 


K 5 


= d ace d bde 


K 8 


^ ace jbde 




K 6 


= d ade d bde 


K 9 


^jo-de jbce 



(20) 
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For the e-scalar interactions a possible basis for N > 4 is therefore 

#i = \Ki 

H 2 = l(K 2 + K 3 ) 

Hi = \K A 

H A = l(K 5 +K 6 ). (22) 
Note that the absence of a d — f type term from the basis follows from the identity 

K s + K 9 = -f^d cde . (23) 

Let us introduce the couplings 

a s = — , a e = and u r = — - , (24) 

47T 4"/T 47T 

and define the corresponding (3 functions for the u T couplings 

^=^A~- (25) 



If we write (with the normalisation of Eq. (112D ) 

4 

\ r H? bcd ~>J2 A7[U r H r (26) 
r=l 

then the /3-functions for the u r couplings are given at one loop bj{§ 

P = Sul + 4N 2 Uxu 2 + 12«2 



A 2 

4(j\T 2 — 4) f 1 
+ — < wim 3 + U1U4 + 2u 2 -u 4 + —(1x3 + 2-U3W4 + 3-u 2 ) 

A, 2 = 12uiw 2 + (2iV 2 + 6)m 2 



— — <! 2m 2 m 3 + 2u 2 u 4 + — (u\ + 6m 3 m 4 + 3m 2 ) 



2(iV 2 -4) 
iV 

j3 U3 = 12mi« 3 + 4-u 2 w 3 + 16-u 2 « 4 

+ {(3iV 2 - 40)u 2 + 6{N 2 - 12)u 3 u 4 + {7N 2 - 96)w 2 } 
(3 U4 = 12mim 4 + 8-u 2 w 3 + 12-u 2 m 4 

+ ^ {4(iV 2 - 14)u 2 + 4(iV 2 - 18)u 3 u 4 -8n 2 } (27) 

where for the moment we suppress contributions from the gauge coupling a s and e-scalar 
Yukawa coupling a e . 



3 Here and for the rest of this section we suppress a factor of I/871" 2 in every one- loop (3- function. 
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Because of the nature of the bare theory, and to explore more easily the supersym- 
metric case, it is natural to consider alternative bases, for example: 



Hi = H 
H 2 = H2 



1 



l^jace jbde _|_ jade j^ ce ^ 

-Jj^ _ l^jaefjbfgjcghjdhe_^_ jaef jbfg jdgh jche^ ^g) 

We shall also see that it is this kind of basis (avoiding use of the d-tensor) that generalises 
most easily to other groups. 



We have 



so that if we write 



then 



— 42 

if 3 = -jyHi — iV"^ 2 ~^ ^ 3 ~~ ^ A 

H 4 = 2Jf 1 + Jf 2 + y#3, (29) 



u r H r = v r H r (30) 

r=l r=l 



Vi = ui-^(u z + u A ) 
2 6 

V 3 = -U4 

2 4 

V4 = N U3+ N U4 ' ^ 
The /^-functions for the v r couplings are given at one loop by 

p vi = 8v 2 + AN 2 v x v 2 + 12v% - 4:Nv!V 3 + 6N 2 V!V 4 + 8Nv 2 v 3 

+ 8Nv 3 v 4 + 8N 2 v 2 v 4 + 10N 2 vj - UNv^s 

(3 V2 = \2v x v 2 + (27V 2 + 6)vj - 4Nv 2 v 3 + 6N 2 v 2 v 4 

- ANv 3 v 4 + 3N 2 vl - 12Nv 2 a s 

(3 V3 = ANv 2 3 + 12 Vl v 3 - Av 2 v 3 - ANv 2 v 4 + (2N 2 + 8)v 3 v 4 

- 2Nv 2 4 - 12Nv 3 a s 

(3 V4 = (fiV 2 + 16)vl + 12 Vl v 4 + 20v 2 v 4 - 2v\ 

- 2Nv 3 v 4 - 12Nv 4 a s + 6a 2 s (32) 

where we have now included the gauge coupling contribution (note that the a 2 terms 
contribute only to /3 V4 ). 
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Another choice of basis (in fact the one employed in Ref. [1]) is 

E x = 5 ac 5 bd + 5 ad 5 bc + 5 ab 5 cd 



H 2 = \(5 ac 5 bd + S ad 5 bc ) - 5 ab 5 cd 
H 3 = H 3 

H A = H 4 (33) 



so that 



and if we write 



then 



#1 = 2(H 1 + H 2 ) 

H 2 = H 2 -2H 1 (34) 
4 4 

v r H r = w r H r (35) 

r=l r=l 



wi = l(vi + 2v 2 ) 

w 2 = l(-v 1 + v 2 ) 

w 3 = v 3 

W4 = V4. (36) 
In this basis the /^-functions become 

Aoi = |[(H2 + l6N 2 )w 2 + {AN 2 - 8)w 1 w 2 - AN Wl w 3 + 26N 2 Wl w 4 

- (2N 2 - A)wl + ANw 2 w 3 + AN 2 w 2 w 4 + 8N 2 w 2 } - UNw^s 
I3 W2 = 1{-8(N 2 + l)wl + 16(N 2 + l) Wl w 2 - WNw^ - 16N 2 W!W 4 

+ (KW 2 - 62)w 2 - 20Nw 2 w 3 + lON 2 w 2 w A - 12Nw 3 w 4 - 7N 2 w 2 } 

- 12Nw 2 a s 

[3 W . A = I6w\w 3 — 8NW1W4 — 28w 2 w 3 — ANw 2 w A + ANwl 

+ (2N 2 + 8)w 3 Wi - 2Nwj - \2Nw 3 a s 

(3 W4 = 64«;i«;4 - 4w 2 w A - 2w 2 - 2Nw 3 w A + (fiV 2 + lQ)wj 

- 12Nw A a s + 6a 2 . (37) 

For the rest of the paper we will use the f-basis; as already remarked the use of H 3 ^ 
means the generalisation to other groups can be carried out easily. 



3.1.1 The fermion contribution 

The contribution of the fermion loop to the scalar anomalous dimension results in a 
contribution of 

A/3 Ui = 8n f I 2 (R)a eUi (38) 
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to each /^-function in Eq. (|27|) . with corresponding contributions to Eq. (|32|) and Eq. (|37j) . 

In Eq. ( |38l) and subsequently we follow the following convention: our fermion repre- 
sentation consists of nf sets of Dirac fermions or 2nf sets of two-component fermions, 
in irreducible representations with identical Casimirs; and the whole representation must 
of course be anomaly free. We will pay particular attention to the case of an adjoint 
representation with n/ = |, which is supersymmetric, and to the case of nf flavours, that 
is rif sets of fundamental two component fermions with rif sets of anti- fundamental two 
component fermions, which is QCD. For the definition of I 2 (R) and more details on group 
theoretic considerations see Appendix [A] 

The 1PI fermion box diagram makes a contribution to the /3-functions (appropriately 
normalised) of the form 

HiA(3 Vi = -4n f a 2 e [Tr(R a R b R c R d ) + Tr(R a R b R d R c ) + Tr(R a R d R c R b ) 

+ Tr(R a R c R d R b ) - Tr(R a R d R b R c ) - Tr(R a R c R b R d )] . (39) 

For a general representation this is not easily expressed in terms of one of our choice of 
bases. In the special case of an adjoint representation (with nf — |), we find that 

HA(3 Vi = a 2 e (-2NH 3 - AH 4 ) (40) 

so that the complete set of /3-functions for the case of an SU(N) gauge theory with an 
adjoint fermion multiplet is: 

/3 Vl = 8v 2 + 4N 2 V!V 2 + I2vl - ANviv^ + 67VVw 4 + 8Nv 2 v 3 

+ 8Nv 3 v 4 + 8N 2 v 2 v 4 + 10N 2 vl - 12N Vl a s + ANv x a e 

(3 V2 = 12 Vl v 2 + (2N 2 + 6)v 2 - ANv 2 v 3 + 6N 2 v 2 v 4 

- \Nv 3 v± + 3A^ 2 V4 - Y2Nv 2 a s + kNv 2 a e 

V3 = ANvl + 12v lV3 - Av 2 v 3 - ANv 2 v± + (2N 2 + 8)v 3 v 4 

- 2Nvl - l2Nv 3 a s + 4Nv 3 a e - 2Na 2 e 
V4 = (fiV 2 + I6)vj + \2v x v A + 20u 2 u 4 - 2v 2 3 

- 2Nv 3 v 4 - 12Nv 4 a s + Qa 2 s + ANv 4 a e - Aa 2 e . (41) 

If we now set V\ = v 2 = v 4 = and v 3 = a e = a s the theory becomes supersymmetric; 
and substituting these values in Eq. (T4T1) we indeed find (5 V1 = (3 V2 = (3 VA = and 

(restoring the 8ir 2 factor) which is identical to the one-loop gauge /3-function j3 s in the 
supersymmetric case. 

Let us consider now the special case of SU(3). In SU(3), (Hi, H 2 , H 3 ) form a basis; 
however if we set v 4 = in Eq. (j4ip then we nevertheless have 

(3 V4 = -2v 2 3 + 6a 2 - 4a 2 . (43) 



10 



This represents a set of contributions to fi Vl2Z which we can identify by using the identity 

H 4 = + H 2 ) + \H 3 . (44) 

Incorporating these contributions into Eq. ( 14TT) we thus find for in SU(3) 







8vf + 36v 1 v 2 + 12v\ - 12v x v 3 + 24^3 - 3§v x a 




+ 


12v\a e — 3v\ + 9a^ — Qa 2 






12v\v 2 + 24i>2 — 12^2^3 — 36v 2 a s 




+ 


12v 2 a e — + 9a^ — Qa 2 






llv 2 + 12viv 3 — <Av 2 v 3 — 3Qv^a s 




+ 


12v 3 a e + 3a 2 s — 8a 2 . 



(45) 

It is easy to check that this set still reduces correctly in the supersymmetric limit. 

For the special case of £77(2), the basis is two dimensional and we have the identities 

2H 1 -H 2 -H 3 = 

2H 1 + H 2 -H A = 0. (46) 

If we choose the basis (77 1; 77 2 ) then we find 

p vi = 8vl + 16viw 2 + 12v% - 2Av x a s + 8v x a e - lQa 2 e + 12a 2 s 

j3 V2 = 12v x v 2 + 14^2 - 2Av 2 a s + 8v 2 a e + Qa 2 . (47) 

Alternatively we could choose the basis (i7 3 , H4) when we find 

(3 V3 = 8v 3 + 24^3^4 — 24t; 3 a s + 8v 3 a e — Aa 2 

(3 Vi = 38vj - 2v\ - Av 3 v A - 2Av 4 a s + 6a 2 + 8v A a e - Aa 2 e . (48) 

With this basis the supersymmetric limit is again apparent; setting w 4 = and v 3 = a e = 
a s we obtain /3 V4 = and (3 V3 = —12a 2 as expected. 

For the fundamental representation of SU(N) we find 

Tr(R a R b R c R d ) = -L [Kt - K 2 + K 3 ] + ~ [K 4 - K 5 + K 6 ] + l -[K 7 + K 8 + K 9 ] (49) 
4iV o o 

and hence for the case of 2rif sets of fermions in the fundamental representation of SU(N) 
(corresponding to QCD with rif flavours), 



HiA(3 v . = 2n f a 



— (77 1 + 77 2 — 77 4 ) — 77 3 



(50) 
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so that the complete set of /3-functions for this case is: 

(3 Vi = 8vf + AN 2 Vl v 2 + \2v\ - AN Vl v 3 + 6N 2 Vl v± + 8Nv 2 v 3 

+ 8Nv 3 v 4 + 8N 2 v 2 v 4 + 10N 2 vj - \2Nv x a 8 + kn s v x a e + 

(3 V2 = I2v x v 2 + (2N 2 + 6)v 2 - ANv 2 v 3 + 6N 2 v 2 v 4 

- ANv 3 v A + 3N 2 vl - 12Nv 2 a s + An f v 2 a e + ^a 2 . 

P V3 = 4Nv 2 3 + I2v x v 3 - Av 2 v 3 - ANv 2 v 4 + {2N 2 + 8)v 3 v 4 

- 2Nvj - \2Nv 3 a s + 4n f v 3 a e - 2n f a 2 e 
(3 V4 = (fiV 2 + 16)v 2 4 + !2v lVi + 20v 2 v 4 - 2v\ 

- 2Nv 3 v 4 -\2Nv 4 a s + §a 2 s + AnfV 4 a e - ^a\. (51) 

It is straightforward to incorporate the fermion contributions in our other choices of bases 
involving Ui or u>j. 

Turning again to the special case of SU(3), and setting v 4 = in Eq. (151]) we have 

p vi = -2v 2 3 + 6a 2 s -^a 2 e (52) 

and incorporating these contributions into Eq. (|5~T|) we thus find for SU(3): 
j3 vi = 8vf + 36v x v 2 + 12t>2 — 12v x v 3 + 24v 2 v 3 — 36v x a s 

Ti f 

+ 4n f v x a e - 3vl + 9a 2 - 2^-a\ 

3 

(3 V2 = 12v x v 2 + 24^2 - I2v 2 v 3 - 36v 2 a s 

Ti f 

+ 4n f v 2 a e - 3vj + 9a 2 - 2^-a\ 

3 

f3 V3 = lltig + I2v x v 3 - 4v 2 v 3 - 3Qv 3 a s 

+ 4n f v 3 a e + 3a 2 - 8^-a 2 e . (53) 

3 

The special case of SU(2) in the fundamental fermion case we leave as an exercise for the 
reader. 



3.2 The general case 

In this subsection we give the results for /3 Vi for a general gauge group. The various 
group invariants are defined in Appendix [A], where results for them for the fundamental 
representations of SU(N), SO(N) and Sp(N) also appear. 

We have derived these results both by substituting in the general expressions that 
follow and by direct calculations with each class of group in the manner described in the 
previous section. 
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We find 



5C 2 A D 2 (RA) + (C A -6C R )D 2 (A)I 2 (R) 2 



^ = " 32n ' 2W\N A -12D 2 {A)(2 + N A ) ^~ l2C ^ 

2 2 192D 2 M) - 80C\N A 
+ SI 2 (R)n fVl a e + S V l + 12 V l- ^(^-3) ^ 

4 f ^-12D 2 (A) + 5C A N A 



27 N A I N A — 3 

72£ 2 (A) 2 - 90C 2 P 3 (A)iV A + 25ClD 2 (A)iV4 ■ 



2hC\N A - \2D 2 (A)(2 + N A ) 
+ Vl [4(l + N A )v 2 - AC A v 3 + QC 2 A v A ] + v 2 (8C A v 3 + 8C> 4 ) 
, _„ 9 5C 2 D 2 (iL4) + (C A - QC R )D 2 (A)I 2 (R) 2 



25C A N A -12D 2 (A)(2 + N 

- \2C A v 2 a s + 8I 2 (R)n f v 2 a e + \2v x v 2 + 2(4 + N A )v 2 
96D 2 (A) - AOC\N A 2 f 12D 2 (A) - 5C\N A 

9C A N A (N A -3) V3VA 27N A \ (N a - 3) 

72D 2 (A) 2 - 90C 2 £> 3 (A)A^ + 25C A D 2 (A)N A } 2 
25C1AT A -12D 2 (A)(2 + 7V A ) r 4 

+ u 2 (-4CUu 3 + 6C> 4 ) 

^ ^ - 25C A N A - UD 2 (A)(2 + ^) { 35 ^ 2 ^ ~ IKfrWW 
+ 4C A £> 2 (iL4)(2 + iV A ) - 16L> 2 (A) 7 2 (i?) (2 + AT A )}a 2 

- 12C A v 3 a s + 8I 2 (R)n f v 3 a e + 12v ± v 3 + AC A vj 

A8D 2 (A)(-1 + N A ) + C A N A (-Q1 + 7N A ) 

+ 2 9ClN A (N A -3) ^ 4 

4 

27C A (N A - 3)N A [25C\N A - 12£> 2 (A)(2 + A^j] 
{l44 J D 2 (A) 2 (2 + AT A )(l + 2iV A ) 
+ 12C1A (A)iV^ [-191 + (-56 + N A )N A ] 

+ Ci^A[-216 J D 3 (^)(-3 + N A )(2 + N A ) + 25(7^(23 + W A )\}vl 

- v 2 (Av 3 + AC A v A ) 

R R 2 , o 5C 2 (C A - 6C R )I 2 (R)N A + 12D 2 (iL4)(2 + iV A ) 2 

A* - ba. + sn; 25Cj7V A - 12D 2 (A)(2 + 7V A ) ^ 

- 2u| — 12C A v±a s + 8I 2 (R)nfV±a e + \2v\V± + 20v 2 v 4 — 2C A w 3 V4 
1152U 3 (A)(2 + iV A ) - 5C 2 [125CjiV A + 4D 2 (A)(98 + N A )\ 2 

6[25C 4 A^-12D 2 (A)(2 + iV A )] ^ 1 J 

The forms taken by C A: r, h(R) and the various invariants D 2 (A) etc for SU(N), SO(N) 
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and Sp(N) are given in Tables [2]|4] in the Appendix. Using Table [2j it is easy to show 
that the results in Eq. (1541) reduce to the results in Eq. (15T|) for the case of SU(N). 



4 The general four-loop results 

The renormalisation constants for the various couplings are defined through 

9s = ^ e Z s g s , 9° e = V e Z e g e , = ^Z Vr ^f¥ r , 



where and T eeee are the one-particle irreducible e-scalar-fermion and four-e-scalar 
Green functions, respectively, the superscript "0" denotes bare quantities, and /i is the 
renormalisation scale. The renormalisation constants associated with the various cou- 
plings satisfy the following relations 

„ _ ^ _ Z** ^ = y/Z* (56) 



ZWy/Z WW ' Z^s/Z^ Z' 

with renormalisation constants as defined in Eq. (fTTj) and Eq. (IT21) 



££ 
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Let us define the f3 functions for the corresponding couplings in the DR scheme: 



a/i z 7r 



e— + 2^L . 
vr Z DR \ da 



77 Z — ' 



DR 

ijklmn 



i,j,k,l,m,n 



77 



8v r 



1 + 2 



dZ? R \ 



Z DR 



7T 



77 



7T 



/3 e (a s ,a e , {w r }) = /U 



d a P 



77 Z, 



d/i 2 77 



5a , 



' Z" fi da e 



77 



-1 



(57) 



o; 



— 6- 



Pijklmn \ ^ J y n 



i,j,k,l,m,n 



77 



77 



77 



77 



(5f 



d/r 7r 



7T A\r 



1 + 2 



2\„ div 



-l 



r'^r 



dv r ' 



2^ p*> 



klmn 



a. 



« e \ J /Vi\ k fV2^ 1 



n 



(59) 

\ 7T / \ 7i / x 77 / \ 77 / \ 77 / \ 77 / 

Here and in the following we do not explicitly display the dependence on the renormal- 
isation scale fx, i.e., a s = a s (fi) etc. Note that in the second line of Eq. ( 1571) . the (9(e) 
terms of j3 e and (3 Vr contribute to the finite part of /3 BR , and similarly for Eqs. (1581) and 
( T59|) . As we will see below, in order to compute the four-loop term of /3 DR one needs (3 e 
to two loops and (3 Vr (r = 1, • • • 4) to one loop. 

For the cases when the fermion representation allows a mass term we introduce the 
fermion mass anomalous dimension, which is defined through 



7 



DR 



/I 2 d 



m DR d^ 2 



-m 



DR 



_7r As — 71 Pe — 77 2 ,Pv r - 



da. 



a. 



77 



da e 



77 



77 



v 3 



77 



C?1V 
77 



(60) 



From this equation one can see that for the four-loop term of 7^ R , the beta functions /3, 
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and (3 Vr are needed to three loops and one loop, respectively, since the dependence of Z® K 
on a e (v r ) starts at one loop (three loops) [1]. The general result for the two-loop order 
f3 e is known [1] and for the three-loop order contributions we find (generalising the QCD 
results from Ref. [2], using a similar calculational setup to the one applied in [1,2], which 
relies on the computer programs QGRAF [15], q2e, exp [16,17] and MINCER [18]): 



A' 



e 
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P 
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P 
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020200 



021010 



P. 



'210010 



P\ 



e 

110101 



211000 



P 
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Pi 



e 
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l-[-QI 2 (R)N A n f - 12N A C R + QN A C A - 10C R + 35C A + 15I 2 (R)n f ] 

ZOO 



63 



1024 



32 
1 



256 



[C\ + QI 2 (R)C R n f - 10C A C R + 2C A I 2 (R)n f + 1QC 2 R ] 
[25N A C A - 3I 2 (R)N A n f - 30N A C R - 5C A - 33I 2 (R)n f + 118C fl ] 



-1-1-27 1 2 (R)n f + 53C A - 30C R ]C A 
IZo 



[hC\l 2 (R)N A + 128D 2 (RA)} 



1024N A I 2 (R) 
17 



= C 

64 • 



15 



512 

— (N A - 1) 

64 V A ' 

^L[3N A C A + lhI 2 {R)N A n f + 6N A C R - 9C A + 3I 2 (R)n f - 50C R ] 



Pi 



010003 



P 



020110 



P: 



e 

110200 



P 



010300 



P\ 
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Pi 
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010021 



e 

110020 



■(96D 2 (A)-7C 4 A N A )C 2 A 



2A576N A 
-^-[35C A + 27 1 2 (R)n f - 18C R ]C A 

IZo 
1 



64 



(N A - 9)C A 

(3N A + 2)(N A -1) 



256 



512 



9 



2048^ 
33 



{32D 2 {A) + 7C\N A ) 



128 



16 



93 



2 



A)ii020 — 5 1 2 ^ 

^020020 = ^[71C A -A2C R -SU 2 (R)n f ]C 2 A 

^010012 = "4096^ 

ne 81 ^3 

Moiom — 512 

A loon = 128 

Monon — 512 

£011200 = ^(Na-1)(N a -2) 

Pe ° 20011 = 153QN a I 2 (R) { ~ 3MCaD2{RA) + 768D ^ A W R ) 

+ N A I 2 (R)[nC\-A2C R C\-m 2 {R)n f C\]} 

£121000 = -7^( 11c a.- &CaC r + 8C r ) 
3 

^010210 = 256^ Na ~ 2 "> Ca 

£110110 = ~~ 32^ 

£011110 = -^(-lO + iV^CU 
9 

£012010 = — 64 ^ 

Mmoio — 32 

£120010 = -^(14^ + 5^(71 

/Wo = -^(^-2^-9/2(^)^1 

£031000 = -^[-l2C A C R + QI 2 (R)C R n f -7C A I 2 (R)n f + 2C 2 A + lQC R ] 

/3 2 e ioioo = -^(-192^ + 47^)^ 

£mioo = ~(5JV A -1)CU 

^10102 = ^i^( 96 ^(^) " UC A N A ) 



17 



•e 

013000 



110002 



A 



112000 



Aoil002 



A 



•e 

010201 



Ao 



011101 



A 



020101 



31 



210001 



Aoi2001 



A 



•e 

021001 



120100 



•e 

030001 



A 



020002 



A 



120001 



+ 



+ 



+ 
+ 
+ 



+ 
+ 



1 



1536iV, 

— (3N A - 
64 1 A 

1 



[-UC A N A + 192D 2 (A)]C A 
5)C A 

(96D 2 (A) - A3C A N A ) 



2048N A 
3 -(5N A -6)C 2 A 



512 



9 
256 



^N A C\ 



~ 768N a I 2 (R) {NaI2{R)[185C3a ~ 207C A^ R ) n f - ^C 2 aCr] 
9Q0D 2 (RA)} 

~ 20A8N A I 2 (R) 1 5C aHR)Na ~ 384C A D 2 (RA) + 96D 2 (A)I 2 (R)} 

-^C 2 2 
128 Ca 

- 7Q8nIi 2 (R) {NaI2{R)[90CaCr ~ U3 ° A ~ ^C 2 A UR)n f ] 
192D 2 (RA)} 
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(3C 2 A - 8C 2 R + \3C A C R ) 



r: ~ 128N a I 2 (R) [{Ca ~ 2C n)i-^D 2 (RA) + C\l 2 {R)N A ] 
- 9I 2 (R)[-8D 2 (RA)+C 3 A I 2 (R)N A }n f } 



~ 18432N A I 2 {R)[-25C\N A + 12D 2 {A){2 + N A )] 
{-3981312D 3 (A) D 2 (RA) + 1132800C A D 2 (RA)N A 
829U0C 3 A D 3 (A)I 2 (R)N A + 12C A D 2 (A)I 2 (R) (18058 - 11N A )N A 
1775C 9 A I 2 {R)N 2 A + 1152C Aj D 2 (A) 2 / 2 ( j R)(586 + 5N A ) 
179712D 3 (RAA)[-25C A N A + 12D 2 (A)(2 + N A )\ 
921QC 2 A [U0C R D 3 (A)I 2 {R)N A + D 2 (A) D 2 (R A) (362 + N A )\ 
75C 8 A I 2 (R)N 2 A [UC R + 27I 2 (R)n f ] 



36C A D 2 (A)I 2 (R)N A [2C R (-16586 + 7N A ) 
9I 2 (R)(206 + 3N A )n f ] - 1990656 D 3 ( A) D 2 (RA)N A 
3456D 2 {A) 2 I 2 {R)[2C R {602 + 13N A ) + 9/ 2 (i?)(2 + N A )n f ]} 

-^ A{ UC A + 5C R ) + ^^(5C - 4C R )D 2 (RA) 
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040000 



A e 3oooo = - UN l Ah{R) (^D 2 (RA) + I 2 (R)N A {11C 3 A - 2A2C 2 A C R 

+ QA0C A C 2 R - 41QC 3 R - 28C 2 A I 2 (R)n f + U4C A C R I 2 (R)n f 

- 10AC 2 R I 2 (R)n f - YlC A I 2 (Rfn) 

+ 48(C A - 2C R ){C A - C R )[2C R -C A + I 2 {R)n f }^}) 

£220000 = ~ 1536 ^/ 2 (fl) [ 3456Z ^K + UR)N A i-^C A - 6A2C 2 A C R 

- 2U8C A C 2 R + 3336C| + 3(247C 2 A + 89QC A C R - 1180C 2 R ) I 2 (R)n f 
+ 24(C A - lQC R )I 2 (R) 2 n 2 f ] - 384D 2 (RA) - 288^24D 2 (iL4) 

+ 2AD 2 {R)n f + I 2 {R)N A {-22C\ + 6C 2 R [6C R - h(R)n f } 

+ 3C A C R [-32C R + h(R)n f ] + C 2 A [81C R + 2I 2 (R)n f ] })c 3 

P™ 000 = ~ 1382AN A I 2 (R) (~ 15552 ^(^) + h{R)N A {\37hhC\ 

- 4C^[13819C i? + 1389I 2 (R)n f ] 

- 8C fl [3483C| - 280I 2 (R) 2 n 2 f + 108C R I 2 (R)n f (-23 + 24C 3 )] 

+ 4C A [12339C£ + 120I 2 (R) 2 n 2 f + 4C R I 2 (R)n f (157 + 1296C 3 )]}) 



192N A I 2 {R)[12D 2 {A){2 + N A ) - 25C\N A ] 
8M<dC A D 2 {RA)I 2 {R)N A n f - 518A0C 2 A C R D 2 (RA)I 2 (R)N A n f 

+ 10368D 2 (RA) 2 (2 + N A )n f + imC 7 A I 2 {R)N 2 A {-A + lllCs) 
150C A N A [C R I 2 (R) 3 N A n 2 f + 2C 2 R I 2 (R) 2 N A n f (23 - 6( 3 ) 

+ 2AD 2 (R)n f (-7 + 2( 3 ) + 8C R I 2 (R)N A (7 + 9C 3 ) 

+ 16D 2 (RA)(-1 + 12Cs)] + 75C A I 2 {R)N 2 A {7I 2 {R) 2 n) 

- AC R I 2 {R)n f {-3\ + 9C 3 ) + 8C 2 R (17 + 6OC3)] 

75C A I 2 (R)N 2 A [I 2 (R)n f (A7 - I6C3) + 4C R (8 + 117Cs)] 

+ 12D 2 {A) (96D 2 (RA)(2 + N A )(-1 + 12(g) 

AC\l 2 {R)N A {2 + 7V A )(-4 + lllC 3 ) 
+ 3C 2 A I 2 (R)N A {4C R (2 + N A )(8 + ll7( 3 ) 
+ I 2 (R)n f [118 + 47N A -16(2 + N A )( 3 ]} 

3C A I 2 (R)N A {7I 2 (R) 2 (2 + N A )n) + 8C 2 R (2 + N A )(17 + 6OC3) 
+ 4C /? / 2 ( J R)n / [134 + 31iV A -9(2 + iV A )C3]} 
+ QC R I 2 (R) 3 N A (2 + iV A )nJ + UAD 2 (R)(2 + N A )n f (-7 + 2C 3 ) 
+ 48C 3 R I 2 (R)N A (2 + N A )(7 + 9( 3 ) 
+ 12C2/ 2 (i?) 2 A^ A n / [262 + 23N A - 6(2 + N A )( 3 ] 
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(Here and elsewhere we denote ((n) by ( n .) We computed the four-loop DRED quantities 
from their DREG counterparts using the indirect method discussed in Refs. [1,19]. It is 
based on the following formulae: 



DK 



7 DR 

Im 



da¥ s 



8v r 



_ ainm DR ^MS^DR ^ dm DR 



d\nm 



MS 



m 



DR 



da™ s 



m 



7T 



!3 Vr dm DR 



m DR dv r 



(62) 



Let us briefly discuss the order in perturbation theory up to which the individual building 
blocks are needed. Of course, the MS quantities are needed to four-loop order; they can be 
found in Refs. [20-23]. The dependence of a s and m DR on a e starts at two- and one- loop 
order [1], respectively. Thus, (3 e is needed up to the three-loop level (cf. Eq. (1621) ). On the 
other hand, both a s and m DR depend on v r starting from three loops and consequently 
only the one-loop term of (3 Vr enters in Eq. ( 1621) . It is given in Eq. (lo^|) . 

For the four-loop analysis we also require the three-loop relations between a s and a^ s 
and between m DR and m MS . The two-loop results were presented in Ref. [1], and the 
three-loop results for the special case of QCD in Ref. [2]. Parametrising the three- loop 
terms by 8^ and 8$, we have 



O'.s 



a 



MS 



1 + 



a 



MS 



7T 12 



C A + 




a 



MS 



a e 1 



7T 7T 



-C R I 2 (R)n f + 5® + . . . 



m 



DR 



■m 



MS 



Ot e 1^ 
7T 4 



a 



MS 



TV 



11 

192' 



C A C 



1 



TV 71 32 



C R (3C A + 8C R ) 



+ (^) 2 ^C R + h(R)n f ]+5^ + ... 



(63) 



where the dots denote higher orders in a^ s , « e , and v r . We find 



7T 



= ^-af s a%(R)n f [2Cl - 3C A C R + 2C 2 R - C A I 2 (R)n f ] 
yb 

+ 7C R I 2 (R)n f ) - T ^(a^) 2 a e J 2 (/2)n / (5Cl + 60C A C R + 6C 2 ) 



+ 



1 



af S (-168C> 4 ^2 - 72C> 4 vi + l2v 3 v 4 C% - 48v 2 v 3 C' 2 A 



9216 



+ 48wit» 3 Cl - A8C A vxv 2 - A8C A N A v t v 2 + 36(7>| + C>| - 72C A u| 
- 2AC A N A v 2 2 - 24C> 2 ) - ^/fvlC A D 2 {A) + ^(aJ^VM^) 



1 



4608 

(«P) 3 
10368 v s ; 



-6C> 3 + 84C> 2 + 36C>! - t> 4 <7 A ) 
[3049<7 A - 416<7 A I 2 ( J R)n / - mC^I^n/] 
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(64) 



a 3 e C R [-10C 2 A + UC A C R + 27C 2 R - 7C A I 2 (R)n } 



384 

+ 39C R I 2 (R)n f - m 2 {Rfn) + 12C% - 36C A C R ( 3 + 24C 2 R ( 3 ] 

- a 2 e C R (jLfiCnV! + 12C R v 2 - 2C A v 2 - C AVl ] 

- 3I 2 (R)C A n f - 19I 2 (R)C R n f - 1Q5C A C R + 144C% 

- 48I 2 (R)C A n f ( 3 + 48I 2 (R)C R n f ( 3 + 72C 2 A ( 3 - 216C A C R ( 3 ]) 

+ a -°R (tt^ [200t; 2 2 + 88N A v 2 2 + 56vf + 167V A ^ 
v 12200 

+ 112N A v 2Vl - C\v\ - l2C\v 3 v A + 176^1 + A8v 3 C A v 2 

- 36C>f + A88C 2 A v A v 2 + 232C> 4 Vi - 48CUuiU 3 ] 

+ ^(af) 2 [2880C% - 168C7 yl / 2 (i2)n / - IUAC A C R - 52C 2 R 



128I 2 (R)C R n f + 1440C% - 4320C A C R ( 3 - 79C\ 
+ ^^(ap) 3 Ci?^A[4354C A + 135C,, + 304/ 2 (i?)n/] 
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4.1 The (3 function and anomalous dimension 

Inserting Eqs. (164")) and (|65|) into Eq. (1621 . we obtain 
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are the one, two, three and four-loop gauge /^-function coefficients calculated in DREG. 
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For the fermion mass anomalous dimension, we find 

74ooooo = 73 + ^ C R C a - §i2 CaC R ~ j^h{R)C A C R n f C 3 

so 99 3 

+ —I 2 (R)C A C 2 R n f + —I 2 ( R yc A C R n} + -I 2 (R)C 2 A C R n f ( 3 

: C R C A l 2 (R)n f - ^^Tl^A^R 



1296 A v ' J 82944 

7X02 = 245 y 6 ^ [2784C A D 2 (A) + 1632C R £> 2 (A) - 53C%JV A 

- \\C R C A N A \C R 

7S11 = -^(53^ + llC R )CiC R 

7120010 = 2tt(3lCA + 18C R )C A C R 

^ = 2048^/ 2 (^) [3 ^ /2( ^^ + 64D2( ^ )] ^ 

7S01 = ^(53^ + 79(7^(71^ 

7oIiio = -^[l7C A + 17I 2 (R)n f -QC R ]C A C R 

7 ^ 011 = ~ 3072A^J 2 (i?) {384jP2(A)/2(jR) ~ 144C aD 2 (RA) 

+ 7V A / 2 ( J R)[-51Ci/ 2 ( J R)n / - 30C*C* + 37C\]}C R 

7 « = - 165888/ 2 (^)iV A ^( 46656D2 (^ 

+ J 2 (i2) iV A {-26505C^ + Cl[355107C R + 23544/ 2 (i?)n / ] 

+ 2C A [11916C| - 65508C ' R I 2 (R)n f - 37 A4I 2 (R) 2 n}} 

+ 12C R [79Q5C 2 R - 7212C R I 2 (R)n f - 224I 2 (R) 2 n 2 f ] 

+ 2592(C A -C R )[5(-11C A -9C R )(C A -2C R ) 

+ 4(5CU - lQC R )I 2 (R)n f }(S 
— 3 

7noo2o = ~2M8^ 3Ca + 11C r) C a c R 

DR _ 7 ^ 2 

72iiooo — 1024 

7 — = ~ 1024iv! 1 J 2 (i?) {ArA/2(jR)h51 ^ /2(jR)n/ ~ 30C2aCr 

+ 37C 3 A ] + 16D 2 (RA)}C R 

- 277Cl/ 2 ( J R)n / ] + 720£> 2 (iL4)}C fl 
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7no2oo = ^(19C R N A + 21C A N A + 95C A + A9C R )C R 

tSioo = -^[^I 2 (R)N A n f + 10C R N A + 3C A N A + 35I 2 (K)n f -U2C R 

- 5C A ]C R 

7oSbio = -^[-51I 2 (R)n f -A2C R + 79C A ]C A C R 

7Soi = ^(U7C R + 2llC A )C 2 A C R 

7?2oioo = -^( c a + qc r+ wc aC r )C r 

.,DR _ 1 ^4^1 

7200011 — ~~ 1024 

.,dr _ J_n?> n 

7200101 — 512 

7S03 = [-41472D 3 (A) + A80C 2 A D 2 (A) + 37C A N A ]C R 

7 ^ 001 = 15MN A I 2 (R) Wa( fl )[~ 63C A ~ 1130^)*/ + 34ClC R ] 

+ 240D 2 (iL4)}C7 fl 
— 1 24 

7P20001 = ^cy<3(3ic A + iscy ^ 

.,DR _ 29 2 

7010120 — 3072 A i? 

7S00 = -^(96^ + 7^)^^ 

.,dr _ 97 3 

7oiom — _ 3072 

.,DR _ *- ^3 r < 

7300010 — 1024 

7300100 = — 512^^ 

DR _ 3 2 

7301000 - ~5i2 A R 

7 ^ 002 = ~ 36864iV A [3360jP2(A) + 1033C ^ Cr 

7S12 = ^^[S8AD 2 (A) + 37C A N A ]C R C A 

7 ™ = -^(^1 + 5^ + 86)^ 

7olooi = -7^( 3iV ^ + 74 )^^ 

_,DR _ 49 

7oiion - 1024 
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Tmoio = ~{^CA + nC R )C A C R 

DR _ ^ nZn 

7200020 ~~ inOA A R 



1024 

-,dr _ 1 n^n 
7201010 — _ 256 



7 2 D oooo2 = ~ 12288AT4 ^ 96D ^ + C A N A \C R C A N A 

7202000 = 5x2^^ 
7201001 — A 

7 ™ = 184 3 2 /^)iv/ ^- 288D2 (^( 1 + 72C ^ 

+ 2C^/ 2 ( J R)jV yl (-1295 + IOO8OC3) 

+ 4 [3^/2(^)^(1544 - 5760Cs) + 2592 J D 2 ( J R)n / (l - 2( 3 ) 

+ 64^/ 2 (i?) 3 iV A ^(-10 + 3Cs) - 2C 2 R I 2 (R) 2 N A n f (571 + IOO8C3)] 

+ C 2 A I 2 (R)N A [2354I 2 (R)n f - 2784I 2 (R)n f ( 3 - 12^(721 + 7704C 3 )] 

+ 2C A / 2 ( J R)iV A [24/ 2 (i?) 2 ^(17-16C3) 

+ C i? / 2 (i?)n / (7444 + 5856C 3 ) + C^(9428 + 71136C 3 )]} 

0/ dr _ }^}_r 2 r 

7ono2o — 3072 

7o ™ = m^ mD2{A) + 37C * ANA]Cn 

7200200 = ^( n a + 3)C a C r 

tSuoo = ^(N A +1)C A C R 

tSoooi = 7f[^-^D 2 {A) + C A N A ]C R 

dr _ 1 r< 2 r 1 

7200110 — 256 

7mioo = ^:(37C A N A + loC R N A + 19C R + 21C A )C R 
oil 

.,DR _ 37 n3 

7oioo30 — 20 48 

7of 2 oo = ^[-^5I 2 (R)N A n f + 39C A N A -38C R N A + 302C R 

- 121I 2 (R)n f - 23C A ]C R 

7 ~ = ^(-^C 2 R -22C A C R + 27C 2 A )C R 
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7040000 



1536I 2 (R)N A [-25C\N A + 12D 2 (A)(2 + N A )\ 



C 



R 



- 16320C 3 A D 2 {RA) I 2 (R) N A n f 

+ 97920C 2 A C R D 2 (RA)I 2 (R)N A n f - 1958AD 2 (RA) 2 (2 + N A )n f 

+ \mC A I 2 {R)N 2 A [-288C R I 2 {R)n f - 15I 2 (R) 2 n 2 f 

+ 2C' (-126 - 432C 3 )] - l00C 7 A I 2 {R)N 2 A {-23 + 240C 3 ) 

+ 100C%I 2 (R)N%{15C R + 7M 2 (R)n f - Q[-132C R + I 2 {R)n f }( 3 } 

+ 25C A N A {13UC 2 R I 2 (R) 2 N A n f + 180C R I 2 (R) 3 N A n 2 f 

- 12n f [5I 2 (R) 4 N A n 2 + 2AD 2 {R){7 - 2( 3 )] 

+ 192D 2 (RA)(-1 + 12C 3 ) + 3C R I 2 (R)N A (368 + 384C 3 )} 

- 12D 2 (A) (-4032D 2 (R)n f + N A {1104C R I 2 (R)(2 + N A ) 

- Q[33QD 2 (R) - 1QC 2 R I 2 (R) 2 (79 + 14^)]% 

+ 180C R I 2 (R) 3 (2 + N A )n) - 60/ 2 ( J R) 4 (2 + N A )n)} 

+ AC A I 2 {R)N A [-2AC R I 2 {R){Al + 12N A )n f - 15/ 2 ( J R) 2 (2 + N A )n) 

+ 2C 2 R (2 + AT A )(-126 - 432C 3 )] + AC\I 2 {R)N A {2 + iV A )(23 - 240C 3 ) 

- 96(2 + N A )[-12C R I 2 (R)N A - 6D 2 (R)n f }( 3 
+ 192D 2 (RA)(2 + N A )(-1 + 12( 3 ) 

+ 4C 2 A I 2 {R)N A {C R {2 + N A )(15 + 792( 3 ) 

+ I 2 (R)n f [184 + 75N A -6(2 + N A )( 3 }}) 

7S01 = -^(23N A + 218)C 2 A C R 

3 

7onmo = ~~^@N A - 2)C A C R 

7mooo = -^(8C A N A -2C R N A -37C A -lhC R )C R 

7X10 = ^(53C A + 11C R )C A C R 

7S00 = - 307 2i^r)n a Cr (~ 57QD2{ra) + ^ R ) N ^ 352C a + QiQC 2 A c R 

- 60A0C A C 2 R + 6992C| - W0C 2 A I 2 (R)n f - 1388C A C R I 2 (R)n f 
+ 3920C£l 2 (i?)n / + lQ8C A I 2 (R) 2 n 2 - 88C R I 2 (R) 2 n 2 

+ A8{C A -C R )[3{8C A -13C R ){C A -2C R ) 

- 2{IQC A - 32C R )I 2 {R)n f + 8I 2 {R) 2 n 2 f ]^]) 

7oSooo = —^[8C A N A -lQC R N A -AI 2 {R)N A n f + A3I 2 {R)n f 

- 50C R + 57C A ]C R 
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,u ™ 3686AI 2 (R)N A [-25C\N A + 12D 2 (A)(2 + N A )} 

C R (-1880064D 3 (A)D 2 (RA) + 536400C A D 2 (RA)N A 

- 391680C 3 A D 3 (A)I 2 (R)N A + AC 5 A D 2 (A)I 2 (R) (22778 - 1UN A )N A 
+ 925^/2(^)^1 + 1152C Aj D 2 (A) 2 / 2 ( j R)(286 + 7N A ) 

+ 8985QD 3 (RAA)[-25C A N A + 12D 2 (A)(2 + JV A )] 

+ 192C^[12240C fl £>3(A)7 2 (i2)^ + D 2 (A) D 2 (RA) (8 198 + 19iV A )] 

- 75C 8 A I 2 (R)N 2 A [10C R + 17I 2 (R)n f ] 

+ 12C A D 2 (A)I 2 (R)N A [-C R (A35A0 - 30N A ) 

+ 3/ 2 (/?)(2634 + 17 N A )n f ] + 3A5Q{-272D 3 (A)D 2 (RA)N A 

+ D 2 (A) 2 I 2 (R)[-2C R (290 + 9N A ) - 13I 2 (R)(2 + N A )n f }}^j 

7£io = ^[3C A -10C R -30I 2 (R)n f ]C 2 A C R 

tEiiio = ^(11^-50)04^ 

7?i2oio = ^Tg(N A + 38)C A C R 

7 ™ = ~ 153QN A I 2 (R) i 384C AD 2 (RA) - 864C R D 2 (RA) 

+ N A I 2 (R)[10C A C R -3C A + 30C A I 2 (R)n f ] 

- 528D 2 (RA)I 2 (R)n f }C R 

7S00 = 3^[ 13 ^ + 8C A I 2 (R)n f + 216C 2 R + 132C R I 2 (R)n f 

- \22C A C R \C R 

7Soo = ^[nC 2 A + 108C 2 R -AlC A I 2 (R)n f -76C A C R + 66C R I 2 (R)n f ]C R 

7 ^ 001 = ~ 4096^/ 2 (fl) ^ 64CaD ^ RA ^ + 32 ^(A)/ 2 (/?) - Cl/ 2 (fl)JV A )C fl 

tSuioi = -^(53^ + 320)^ 

7^200 = -^(5^1 + 21^ + 46)^ 

7?vIioo = -^(^1 + 31^ + 22)^ 

7Soo = ^(N A -W)C R 
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where 



73 



1 

256 



Cj(-ifi- 336C3) + CjC. 4 (if%316C: 
34045 



/ 15349 



^ 2 ^of 34U45 AAn ^\ „ ^,/70055 1418 A .„ nf .\ 

+ CflCl ( gg- - l52Cs + 440C 5 j + C R C\ {-^ + — Cs - 440Cs j 



+ C R I 2 (R)n f 
+ C 2 R C A I 2 (R)n f 
+ C R C 2 A I 2 (R)n f 



280 



+ 552C 3 - 4800 



8819 

65459 
162 



368Cs - 264( 4 + 80Ce 



2684 



C 3 + 264C 4 + 400C E 



+ ^/ 2 ( J R) 2 nj(^-l60C 3 + 96C4 

+ ^(^)(-^^C3) + ^(-3 2 + 24 0C3) 



+ C R C A I 2 (R) 2 n 2 



1342 



+ 160C3-96C4 + 



n f D 2 {R) 



64 - 480C; 



(69) 



5 The four-loop supersymmetric case 

Our conventions are such that substituting in the above equations the results of Tables EH 
corresponds to a gauge theory with rif sets of Dirac fermions transforming according to 
the fundamental representation, or rif sets of fundamental two component fermions with 
rif sets of anti- fundamental two component fermions. 

To extract the supersymmetric case we must make the replacements 



Cr 

h{R) 

n f 
D 2 {R) 

D 2 (RA) 

D 3 (RAA) 

a e 

Vi 



C A 
C A 



D 2 (A) 

D 2 (A) 

D 3 (A) 

v 3 = a s 

v 2 = Vi = 0. 



(70) 



With these substitutions we can compare our results for (3 S with the four-loop results for 
the gauge /^-function /3f YM of SQCD which was presented in Ref. [6]: 



0. 



SYM 

s 



7T 



-c A + 3 -c 2 A ^ + 2 -±c A (^)\ 

4 8 A % 64 A \n ' 



51 
128 



C\l- 

7T 



0(a°) . (71) 



We indeed find that using Eq. (1701 in Eq. (1571) precisely reproduces Eq. (1711) . 
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We have also checked that in the same supersymmetric limit, Eq. (|6T|) reproduces the 
three-loop SQCD /^-function. 

Turning now to the case of softly-broken supersymmetry, there exists an exact result 
for 7 m [8]: 



SYM 

/ m 



da. 



SYM 



a. 



(72) 



SYM 

/ m 



51 
— ( 

32 



-) +^c%[- 



+ 0{a 



(73) 



whence it follows that 

4 Ca + A Ca V + 64 C 
Using Eq. ([70]) in Eq. (I60I) precisely reproduces Eq. (1731 in similar fashion. 

The invariant D 3 (A) does not feature in either calculation, and the dependence on 
D 2 (A), Na, C3; C4 an d Cs & U cancel, although they appear in individual terms. It is 
tempting to speculate that this absence of higher order invariants and transcendental 
numbers (other than n) is related to the existence of the NSVZ scheme, in which the 
gauge /3-function for any simple gauge group is given (in the supersymmetric case without 
matter fields) by the expression [24], [25] 

-1 



P, 



NSVZ 



1 



2?r 



(74) 



which is manifestly free of them to all orders. It is natural to conjecture that the same 
property holds in the DRED scheme. For discussion of the relationship between [3f svz 
and ^F* see Ref. [7]. 



6 Discussion 



In this paper we have applied DRED to gauge theories with gauge groups SU(N), SO(N) 
and Sp(N), and calculated both the gauge /3-function and the mass anomalous dimension 
to the four-loop level. These calculations required careful treatment of the evanescent 
Yukawa and quartic couplings of the e-scalar. In the supersymmetric limit we explicitly 
verified that the /3-function for the evanescent Yukawa coupling reproduces the gauge 
/5-function through three loops. 

The results for /3° R and 7° R in the special case of QCD as described in Ref. [1] and [2] 
are easily obtained from the results of this paper by specialising to the SU (N) case, and 
setting N = 3 and v± = 0, with the fermions in the fundamental representation. 

Predictions based on theories with low energy supersymmetry require careful consider- 
ation of the transition between the MS and DR renormalisation schemes. If, for example, 
the decoupling of supersymmetric particles is carried out in several steps (as in split su- 
persymmetry, for example [26]) then it is essential to take into account the evanescent 
couplings (for a recent discussion and treatment of the running of a s and mj, in the MSSM, 
see Ref. [27]. 
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A Group Theory 



We consider a gauge group Q with generators R a satisfying^] 

[R a ,R b ] =if abc R c . (75) 

We work throughout with a fermion representation consisting of n/ sets of Dirac fermions 
or 2nf sets of two-component fermions, in irreducible representations with identical Casimirs, 
using R a to denote the generators in one such representation. Thus R a R a is proportional 
to the unit matrix: 

R a R a = C R .I (76) 
For the adjoint representation we have 

CA^ab = facdfbcd- (77) 

I 2 (R) is defined by 

Tr[R a R b ] = I 2 {R)8 ah . (78) 

Thus we have 

C R d R = I 2 {R)N A (79) 

where Na is the number of generators and cLr is the dimensionality of the representation 
R. Evidently I 2 {A) = Ca- The fully symmetric tensors d a R bcd and d a A bcd are defined by 

d abcd = ±.T T [R(*Rb R c R d)^ 

d abcd = _^ F {a F b F cpd)^ ^ 



where 
and 



(F a ) bc = %f ac (81) 

R( a R b R c R d) = R a R b R c R d + R a R b R d R c + R a R c R b R d 

+ R a R c R d R b + R a R d R b R c + R a R d R c R b , (82) 

(similarly for F^ a F b F c F d) ). 

The additional tensor invariants occurring in our results for (3 S and 7 m are defined as 

D 2 (A) = d a A bcd d a A bcd 

D 2 (RA) = d abcd d a A bcd 

D 3 {A) = d abcd df f d abef 

D 3 (RAA) = d a R bcd d A ef d a A bef . (83) 

In table [2}|4] we present results for the various tensor invariants for the groups SU(N), 
SO(N) and Sp(N), when the fermion representation R is the fundamental representation. 
In each case the constant b reflects the arbitrariness in the choice of normalisation of the 
generators (see Eq. (fT8|) for SU(N)). 

4 Useful sources for some of the material in this section have included Refs. [11,20,28]. 
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B The groups SO{N) and Sp{N) 



In this section we derive explicit expressions for the /3-functions for the e-scalar quartic 
interactions for the groups SO(N) and Sp(N). These may also be derived from Eqs. fl54l) 
using tables [3l HI 



B.l The case Q = SO{N) 

Let us consider SO(N). The defining representation of the generators of SO(N) is given 
by the set of matrices 

(M[ij]) k i = -i(5ik5ji - 5 u 5 k j), (84) 

satisfying the algebra 

[M {ij] ,M [kl] ] = -i (8 jk M m - 8 ik M m - 8 fi M m + 8 U M M ) (85) 

or 

[M[ij] , Mm] = if[ij][ k i\[mn]M[ mn ] (86) 
where the structure constants are given by 

f[ij][ki][mn] = \ [8ik8j m 8in + ■■■ (7terms)] (87) 

such that they are antisymmetric in ij, kl, and mn exchange. With this normalisation 
of the generators it is straightforward to show that the adjoint quadratic Casimir Ca is 
given by 

f[ij][kl}[mn}f[ij}[kl}[ P q] = CA{8 mp 5 nq — 8 mq 8 np ) (88) 

with 

C A = 2(N-2). (89) 

We will, however, present results for an arbitrary normalisation of the generators such 
that 

C A = b(N-2), (90) 

where b is a constant. Useful checks on our calculations will be provided by the isomor- 
phisms 

SO(3) ^ (91) 

SO® , ™& (92) 

which mean that the Lie algebras of 50(3) and SU(2), and of 50(6) and SU(A) are 
identical. Note that to compare our result for 50(3) with the corresponding result for 
SU(2) (where with the conventional normalisation we have Ca = 2) we will need to set 
b = 2, while to compare 50(6) with SU(4) we will similarly need to set 6=1. 
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The basis for 4-tensors for SO(N) for N > A has 7 = 6 and can be chosen to bd^|(we 
adopt a shorthand notation with [M2] — > i etc.): 

Pi = SijSki, 
P3 = SuSkj, 

P4 fijmfklmi 
P5 fikmfjlmt 

P 6 = tvlFiFjFkFi], 

(93) 

where (Fj) mn = f min . 

Some useful identities for reduction of various 4-tensors to the basis are as follows : 

Tr Wi Pj Pp Fk Pi Pp) = y(iV-4)(-2F 1 + 2P 2 + P 3 ) 

+ ^(iV-8)(P 4 -2P 5 )-^(iV-2)P 6 (94) 

Tr [F{ F m Fj F n ](F k Fi) mn = ^{N - A){-2P 1 + P 2 + P 3 ) 

- ^(AT-8)F 4 + ^(iV-8)F 5 . (95) 
It is interesting that this doesn't involve the basis element P$. 

[Pi Fj F m Pn] Tr [F k F t F m F n ] = j(N - 2)(N - 4)(6P - P 2 - P 3 ) 

b 3 

+ — (N 3 -6N 2 + 16N-24)P 4 
8 

- ^N(N-6)P 5 
b 2 

+ —{N 2 — 6N + 20)P 6 (96) 

Tr [Fj F m Fj F n F k F m Fi F n ] = j ( N - 2) {N - 4) (P x + P 3 ) 

b 3 

+ —(N 2 — IAN + 32)(— 2P 4 + P 5 ) 
8 

- ^(iV-8)P 6 (97) 



5 As will become clear our results will not be applicable to the special case N = 8, which we will not 
consider further 
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The results for the one-loop /3-functions in a basis as in Eq. (1281) and with v i defined 
as for SU(N) follows: 

p vi = 8v 2 + 2(N 2 -N + 2)vtv 2 -Ab(N- 2)v 1 v 3 + 6b 2 (N -2) 2 v l v i 

+ l2v 2 + 8b(N-2)v 2 v 3 + 8b 2 (N-2) 2 v 2 v 4 

+ 206 4 (iV - 2)(JV - A)v 2 4 + 166 3 (iV - 4)v 3 v 4 - 12b(N - 2) Vl a s , 

Pv2 = (N 2 -N + 8)v 2 + l2v 1 v 2 -4b(N-2)v 2 v 3 + Qb 2 (N-2) 2 v 2 v 4 

- 8b 3 (N - A)v 3 v 4 + 6b 4 (N -2)(N- A)vj - 12b(N - 2)v 2 a s , 
(3 V3 = Ab(N - 2)v 2 3 + 12v 4 v 3 - Av 2 v 3 - Ab(N - 2)v 2 v 4 

+ 2b 2 (N 2 -6N + 20)v 3 v 4 - Ab 3 (N - A)v\ - l2b(N - 2)v 3 a s , 
[3 V4 = |6 2 (3A^ 2 - 28^ + 140)^ + 12^^ + 20^4-2^ 

- 2b(N - 2)v 3 v 4 - \2b[N - 2)v 4 a s + 6a 2 s . (98) 

If we set iV = 6 and b = 1 in Eq. ( 1981) . we reproduce precisely the results obtained by 
setting N = 4 in Eq. (1321) . Similarly, if we set iV = 3 and b = 2 in Eq. (1981) . we reproduce 
precisely the results obtained by setting N = 2 in Eq. (|32|) . setting v 3 = v 4 = in both 
cases. 



B.l.l The fermion contribution 

As in the SU(N) case, the fermion loop contribution to the scalar anomalous dimension 
results in a contribution of A/3 Vi = 8rifI 2 (R)via e to each /^-function in Eq. (1981) . The 1PI 
fermion box diagram makes a contribution to the /^-functions (appropriately normalised) 
of the form as Eq. (1391) ; for the adjoint representation this becomes: 

HAf3 Vi = a 2 e (-2b(N - 2)H 3 - 4i7 4 ), (99) 

and hence for the complete /^-functions including a single two-component fermion in the 
adjoint representation we have from Eq. fl98l) : 

P n = 8vf + 2(N 2 -N + 2) Vl v 2 - Ab(N-2)v 1 v 3 + 6b 2 (N- 2) 2 Vl v 4 

+ 12v 2 + 8b(N - 2)v 2 v 3 + 8b 2 (N - 2) 2 v 2 v 4 + lQb 3 (N - A)v 3 v 4 

+ 2%\N - 2)(JV - A)vj - 12b(N - 2) Vl a s + Ab(N - 2) Vl a e , 

(3 V2 = (N 2 - N + 8)v 2 + 12v 4 v 2 - Ab(N - 2)v 2 v 3 + Qb 2 (N - 2) 2 v 2 v 4 

- 8b 3 (N - 4)v 3 v 4 + 6b\N -2)(N - 4)v 2 4 

- 12b(N - 2)v 2 a s + Ab(N - 2)v 2 a e , 

(3 V3 = Ab{N-2)v 2 +12v 1 v 3 -iv 2 v 3 -ib{N-2)v 2 v 4 

+ 2b 2 (N 2 - 6N + 20)v 3 v 4 - Ab 3 (N - A)v\ - l2b(N - 2)v 3 a s 

+ Ab(N - 2)v 3 a e - 2b(N - 2)a 2 e , 

(3 V4 = ±b 2 (3N 2 -28N + U0)vl + l2 Vl v 4 + 20v 2 v 4 -2v 2 3 

- 2b(N - 2)v 3 v 4 - 12b(N - 2)v 4 a s + 6a 2 s 

+ Ab(N - 2)v 4 a e - Aa 2 e , (100) 
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when it is once again easy to extract the supersymmetric result by setting v\ = vi = v± = 
and t>3 = a e = a s . 

For a single two-component fermion in the fundamental representation, (and for N 7^ 
8) we find that 

Tr[MiMjM k Mi] = — i- \-b 2 (P 1 + P 2 + P 3 ) + b(2P 4 - P 5 ) + Pel (101) 

1\ — o L 

and hence we find for fermions in the fundamental representation a contribution to the 
/3-functions (for n/ flavours) of the form 

= [86 2 Cffx + H 2 ) - 2b(N - 10)H 3 - AH 4 ] . (102) 

It is straightforward to incorporate these contributions into Eq. ( |98i) in the same manner. 



B.2 The case Q = Sp{N) 

We will here be considering the unitary symplectic group. The generators of Sp(N) satisfy 

JR a J = {R a ) T (103) 

where 

1 - (-/ 0) < 104 > 

and I is the unit matrix. Evidently iV must be even. For the case iV = 2 it is easy to 
check by explicitly constructing R a to satisfy Eq. H103[) that Sp(2) = SU(2). Another 
useful check on our calculations will be provided by the isomorphism 

SO(5) = ™ ( 10 5) 

^2 

If we write N = 2n, the generators may be written as L a p, where an infinitesimal 
group element S may be written 

S = 1 + i yj a a pL a p (106) 

a/3 

where a a p = ap a and 

L a /3 = L_p_ a , a,/3 = ±l,±2,---±n (107) 

(Thus the correspondence Sp(2) ~ SU{2) is L\\ ~ J3, £1-1, L-1,1 ~ J± = J\ ± ^2-) 
They obey the commutation relations 

[L aj g, L 7( j] = (8pry8pL a s — SasdaL^p + 8pg6pL a y — S a 7y8 a Lp S ) (108) 
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where 



(109) 



and a = —a etc. 
We find: 

P V1 = 8vl + 2(N 2 + N + 2)v x v 2 -2b{N + 2) Vl v 3 + p 2 (N + 2) 2 vm 

+ 12vj + Ab(N + 2)v 2 v 3 + 2b 2 (N + 2) 2 v 2 v A 

+ \b\N + 2)(N + A)v\ + 2b\N + A)v 3 v A - 6b(N + 2) Vl a s , 

P V2 = (N 2 + N + 8)v% + 12v lV2 -2b(N + 2)v 2 v 3 + \b 2 {N + 2) 2 v 2 v i 

- b\N + A)v 3 v 4 + \b\N + 2)(N + A)v 2 A - 6b(N + 2)v 2 a s , 
(5 V3 = 2b(N + 2)vl + Uvm - Av 2 v 3 - 2b(N + 2)v 2 v 4 

+ \b 2 (N 2 + 6N + 20)^4 - \b\N + A)v\ - 6b(N + 2)v 3 a s , 

(3 VA = |6 2 (3iV 2 + 28A^ + 140)t»2 + 12t» 1 t>4 + 20^4-2^ 

- b(N + 2)v 3 Vi-6b(N + 2)v A a s + 6a 2 s . (110) 

Using table H] in Eq. fl54l) leads to the same results. 

Setting N = 2, b = 1 and v 3 = i> 4 = in Eq. (jllOp above we indeed find agreement 
with the corresponding results for SU(2), from Eq. fl32l . 



B.2.1 The fermion contribution 

The fermion loop contribution is similar to the the SO(N) case. From the scalar anoma- 
lous dimension we get a contribution of A(3 V . = 8rifI 2 (R)via e to each ,5-function in 
Eq. (THUD. 

In the case of the adjoint representation, by similar algebra to that leading to Eq. (j4"Ul) . 
we obtain for the 1PI fermion box diagram a contribution: 

HAl3 Vi = a 2 e (-b(N + 2)H 3 - AH 4 ), (111) 



38 



and for the complete /3-functions from Eq. (11 101) : 



(3 V1 = 8vf + 2(N 2 + N + 2)v x v 2 -2b(N + 2) Vl v 3 + \b 2 {N + 2) 2 Vl v 4 
+ 12v 2 + Ab(N + 2)v 2 v 3 + 2b 2 (N + 2) 2 v 2 v A 



+ §£> 4 (N + 2)(JV + A)v\ + 2b 3 (N + A)v 3 v 4 - 6b(N + 2) Vl a s 

+ 2b{N + 2)v 1 a e , 

P V2 = (N 2 + N + 8)v 2 2 + 12 Vl v 2 - 2b(N + 2)v 2 v 3 + \b 2 (N + 2) 2 v 2 v 4 

- b 3 (N + A)v 3 v A + lb\N + 2)(N + A)v 2 A - 6b(N + 2)v 2 a s 
+ 2b(N + 2)v 2 a e , 

P V3 = 2b(N + 2)v 2 3 + 12^3 - Av 2 v 3 - 2b{N + 2)v 2 v 4 

+ \b 2 (N 2 + QN + 20)v 3 v 4 - ±b 3 (N + A)vj - Qb(N + 2)v 3 a s , 

+ 2b(N + 2)v 3 a e - b(N + 2)a 2 e , 

(] V4 = ±b 2 (3N 2 + 28N + U0)vl + l2 Vl v 4 + 20v 2 v 4 -2vl 

- b{N + 2)v 3 v 4 - 6b(N + 2)v 4 a s + 6a 2 s 

+ 2b(N + 2)v 4 a e - Aa 2 e , (112) 

when it is once again easy to extract the supersymmetric result by setting v\ = v 2 = v A = 
and v 3 = a e = a s . Setting = 4 in Eq. (11121) we reproduce precisely the results of setting 
iV = 5 in Eq. (1981) . in accordance with Eq. (11051) ; a good check of our calculation. Also, 
setting N = 2 and v 3 = i>4 = in Eq. (11121) we reproduce the results of setting N = 2 
and v 3 = i> 4 = in Eq. (I3"2"j) . 

For a single two-component fermion in the fundamental representation, we find (again 
reverting to a shorthand single index notation) 

Tr[L a LpL 7 L s ] = [- j(Pi + P 2 + P 3 ) - \{2P A - P 5 ) + P 6 ] , (113) 

(where the P-basis is defined in the same way as for SO(N) in Eq. (1931 ) and hence a 
contribution to the /3-functions (for rif flavours) of the form 

HiAp Vi = -f-§ [2b 2 (Ht + H 2 ) - b(N + 10)H 3 - AH A ] . (114) 

iv + o 

It is straightforward to incorporate these contributions into Eq. (IllOp in the same manner. 

We can also check this result using the identity Sp(2) = SU(2); setting N = 2 and 
b = 1 in Eq. (11141) and using Eq. (T46|) we find agreement with the result of setting N = 2 
and using Eq. (H6l) in Eq. (|50|) . 

Remarkably, the all /3-functions for the Sp(N) case (including (3 S , f3 e ) together with 7 m 



6 Of course we cannot use Eq. (|105[) as a check here because the fundamental representation is different 
for the two groups. 
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can be derived from the corresponding SO(N) versions by a series of simple substitutions: 



b - 


♦ \b 


N - 


-> -N 


a s - 


-> -a s 


a e - 


-> -a e 


v 3 - 


-> -v 3 


n f - 


-> -71/ 



(115) 
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Group 


SU{N) 


c A 


bN 


C R 


b N*-l 
u 2N 


h(A) 


bN 


h(R) 


b 
2 


N A 


N 2 -l 


D 2 (A) 


^(N 2 -1)(N 2 + 3Q)N 2 


D 2 (RA) 


^N(N 2 -l)(N 2 + 6) 


D 2 {R) 


^(N 2 - 1)(18 - QN 2 + N*) 


DM) 


^N 2 (N 2 - 1)(324 + 135iV 2 + iV 4 ) 


D 3 (RAA) 


J^V 3 (iV 2 -l)(51 + Ar 2 ) 



Table 2: SU(N) Group invariants (here R is the fundamental representation). 
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Group 


SO{N) 


Ca 


b(N - 2) 


w H 


^-(N - 1) 


hi A) 


b(N-2) 


h{R) 


b 


N A 


±N(N-1) 


D 2 (A) 


^N(N - 1)(N - 2)(-296 + 138iV - 15iV 2 + N 3 ) 


D 2 {RA) 


^N(N - 1)(N - 2)(22 -7N + N 2 ) 


D 2 (R) 


^N(N-1)(4-N + N 2 ) 


D*(A) 


^(N -2)(N- l)7V(-29440 + 23272iV - 7018iV 2 + 9717V 3 - 47iV 4 + 2N 5 ) 


D 3 (RAA) 


^N(N -2)(N- 1)(2048 - 1582iV + 387iV 2 - 317V 3 + 2iV 4 ) 



Table 3: SO(N) Group invariants (here R is the fundamental representation). 
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Group 


Sp(N) 


Ca 


b(N + 2) 

y 1 ^ J 


w H 


-(N + 1) 


hi A) 


b(N + 2) 


h(R) 


b 
2 


N A 


±N(N + 1) 


D 2 (A) 


^N(N + 1)(N + 2)(296 + 1387V + 15iV 2 + iV 3 ) 

f Do 


D 2 {RA) 


^N(N + 1)(N + 2)(22 + 7N + N 2 ) 


D 2 (R) 


^N(N+1)(4 + N + N 2 ) 


D*(A) 


55296 + 2 )(^ + 1)^(29440 + 23272iV + 7018iV 2 + 9717V 3 + 477V 4 + 27V 5 ) 


D 3 (RAA) 


55296 N ( N + 2 )( N + 1)(2048 + 1582iV + 387iV 2 + 31iV 3 + 27V 4 ) 



Table 4: Sp(N) Group invariants (here R is the fundamental representation). 
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